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Abstract 


In this paper, we study and compute the inverse of matrices with 
parametric entries. We demonstrate that the Gauss-Jordan method 
can be extended to compute the inverse of parametric matrices, offering 
a powerful tool for solving systems of linear equations and analyzing 
parametric systems. Using this new expansion (so-called Gauss-Jordan 
systems) and also utilizing linearly dependency systems for linear sys- 
tems involving parameters [4, 5], we introduce the notion of an inverse 
matrix system for a parametric matrix. In doing so, we decompose 
the space of parameters into a finite partition and for each partition, 
we give the corresponding inverse matrix without applying Grobner 
systems. We also present an algorithm for computing an inverse system 
for a given parametric matrix. All mentioned algorithms have been 
implemented in Maple, and their efficiency and behavior have been 
experimented on a set of benchmark matrices. 
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1 Introduction 


The inverse matrix computation is a fundamental operation in linear alge- 
bra with profound applications across various fields, including engineering, 
physics, computer science, and statistics. The inverse matrix is utilized as 
a powerful tool for solving systems of linear equations, providing insights 
into the underlying structures of mathematical models. Among the myriad 
techniques available for computing the inverse of a matrix, the Gauss-Jordan 
method stands out as a versatile approach that has found widespread use 
due to its ease of implementation. 

Often, matrices with parametric entries, where elements of the matrix 
are expressed as parameters, arise in modeling real-world problems. This 
paper explores the application of the Gauss-Jordan method to compute the 
inverse of matrices with parametric entries. We demonstrate how the Gauss- 
Jordan method can be extended to handle matrices with parametric entries, 
empowering researchers and practitioners to solve complex problems in a 
variety of domains. The following example shows that the common approach 
for computing the inverse of matrices may not be used in parametric cases. 


Example 1 Let 


—b 1 a+l1 
A= 0 Cc b+1 
—-1 3+ce 1 


where a,b and c are parameters in the complex numbers. The matrix in- 
verse of A computed by the function LinearAlgebra:-MatrixInverse(A) 
of Maple is equal to 


bc+3 64+3 ac+3 at+c+2 ac—b+c—1 
b2c+act+3 b24+2b+c—1 b*%ce+ac+3 b?+2b+c—1 b2c+ac+3 b2+2b+c—1 
Aw} a b+1 a—b+1 (b+1)b 
b2c+ac+3 b24+2b+c—-1 b2ce+ac+3b24+2b+c—1 b2c+ac+3 b?4+2b+c—1 
c bc+3 b—-1 be 
b2c+ac+3 b2+2b+c—1 b2c+ac+3 b2+2 b+c—1 b2c+ac+3 b2+2b+c—1 


Considering A = b?c + ac + 3b? + 2b+c¢—1 as the determinant of A the 
matric A! is defined when A #0. For instance if a=1, b=0, andc=0 
then A and A7! are as follow: 
0. 2 
0 01 


Alta=1,b=0,c=0} = 
hs 3 | 
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3 -5 -1 
A™|fa=1,b=0,c=0} = | 1 -2 0 
y 1 | 
However, A~' undefined for the values a = 1,b = } and c = 0, since 


A = bc+ac+ 3b? + 2b+c—1=0 and Al ta=1,b=3,c=0} is the following 
Matrix with rank 2. 


=1/3 1 2 
Al ig=-95=-tye=0) = 0 0 4/3 
-1 3 1 


In fact, one may be interested in identifying all possible parametric relations 
between a, b, and c where A is invertible or singular. With more details, we 
want to partition V(A) and also C?\V(A) into subvarieties. For this purpose, 
we introduce the inverse matrix system for parametric matrices and present 
an algorithm for computing it. In doing so, we can use two ideas; Grobner 
systems and Gauss-Jordan systems based on linearly dependency systems. 


2 Grobner Systems 


The concepts of Grobner systems together with their algorithms were intro- 
duced by Weispfenning [13] in 1992. Grobner systems can be considered as 
an extension of Grobner bases to polynomials with parametric coefficients. 
Roughly speaking, for a parametric ideal, by computing its Grobner system, 
we partition the space of parameters into a finite set of cells, and for each 
cell, we provide a set of polynomials. For any values of parameters, we find 
out the cell containing these values, and the corresponding polynomial set is 
a Grobner basis of the ideal for the given values of parameters. 
Throughout this article, we consider R = K[xj,..., 2] the polynomial 
ring in terms of 21,...,%m over a field K. Let Z = (fi,..., fg) C R be the 
polynomial ideal generated by the f;’s. For any f € R, the leading monomial 
of f, denoted by LM_(f), is the greatest monomial (w.r.t. the monomial 
ordering <) appearing in f and its coefficient is the leading coefficient of f 
which denoted by LC.(f). The leading term of f w.r.t. ~< is the product 
LT.(f) = LC.(f)LM.(f). The leading monomial ideal of T is defined 
to be LM.(Z) = (LM.(f) | f € Z). A finite subset {91,...,gm} C T is 
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called a Grébner basis for T w.r.t. ~ if LM.(Z) = (LM.(g1),..., LM2(gm)).- 
Grobner bases together with the first algorithm to compute them were 
introduced by Buchberger in his Ph.D. thesis under the supervision of 
Grobner [1]. For more details, we refer the interested readers to [2]. 

Using these notations, we recall the definition of Grébner bases for 
ideals containing polynomial generators with parametric coefficients, namely 
Groébner systems. Roughly speaking, Grébner systems can be considered as 
an extension of Grobner bases for polynomial ideals over fields to polynomial 
ideals with parametric coefficients. For this purpose, let us consider S = 
K[a, x] as a polynomial ring with parametric coefficients where x = 71,...,2n 
is a sequence of variables and a = aj,...,@m is a sequence of parameters. 
Also, {x} and {a} are disjoint sets ({x}M {a} = @) which means that x; 4 a; 
for any i = 1,...,n and j = 1,...,m. Thus any monomial in S is in the 
form x® = zf'--- a9" where a = (a,...,Qn) € N”. Additionally, a term 
in S is a’x® which a7 € K{a] is a monomial as a coefficient of x®. Let 
us consider <, and <, as two monomial orderings on the variables and 
parameters, respectively. To define and compute Grobner systems, we also 
need to recall a product order to specify an ordering on S. The product 
ordering <x,q, the product of two monomial orderings <, and ~a, is defined 
as follows: For any y,6 € N™ and a, 8 € N”, we can write a’x® <x a aox® if 
either x® <, x? or 1° = xPand a? <, ao}. Furthermore, an specialization 
of parameters is a morphism o : K[a] + K which in K denotes the algebraic 
closure of K. So, for any m—tuple a = t1,...,tm and for each polynomial 
f € S one may be write o(f) = fla=t,,...tm Where o(a;) = t;. Moreover, a 
parametric polynomial ideal is a type of polynomial ideal defined by a system 
of parametric polynomial equations. This ideal is generated by polynomials 
with parameter coefficients. Additionally, a parametric linear polynomial 
refers to a polynomial where the power of the variables x;’s is at most one, 
and the parametric coefficients can be polynomials of any power in K[a]. 

Now, let’s review the definition of a Grébner system for parametric 
polynomial ideals. The concept of Grdébner system was introduced by 
Weispfenning in [13]. He proved that any parametric polynomial ideal has 
a Grébner system [13, Proposition 3.4 and Theorem 2.7] and described an 
algorithm to compute it [13, Theorem 3.6]. Since then, several different 
algorithms have been designed to compute the Grobner system of an ideal, 
each of which has advantages and faults [4, 5, 3, 6, 7, 8, 9, 10, 11, 12]. 


Definition 1 Let F = {fi,..., fe} C S andG = {(Ni, Wi, Gi), be a 
finite triple set where N;,W; C Kia] and G; C S are finite fori =1,...,2. 
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The set G is called a Grobner system for (F) w.r.t. to Xx over VCR” if 
for any i we have 


e 0(G;) is a Grébner basis of (o(F')) with respect to <x, for any special- 
ization o : K[a] > K satisfying (Ni, Wi); ( if o(p) =0 for all p € N; 
and o(q) £0 for all q € W; then o satisfies (Ni,Wi) C Kla] x K[a]). 


eVC (is V(N;) \ V(W;) which V(N) is the common zeros of a set of 
polynomial generators for N. 


Any triple (Ni, Wi, Gi) is said a branch (segment) of the Groébner system G 
for each 1 < i < &. Moreover, if VY = K” then G is generally called a 
Grébner system of F’. Also, each pair (N;,W;) is called a specification or a 
parametric constraint ( N; and W; are called the null and non-null condition 
sets, respectively). Additionally, a parametric constraint (NV,W) is called 
inconsistent if and only if for any polynomial f € W, f belongs in the radical 
ideal generated by N. An effective technique for deciding the consistency of 
parametric constraints involves using an efficient combination of the ICheck 
and CCheck algorithms (introduced in [8]), along with the incorporation of 
Rabinovitch’s trick for testing radical membership in the worst-case scenario; 
this is achieved by introducing a new variable. The interested reader can be 
referred to [8] for more details. 


Example 2 Let F = {(a — b)ry + (b— 1)y, x? + (2 — b)y} C K[a,b, 2, y] 
where x,y are variables and a,b are parameters. We consider the monomial 
orderings z ~<drt Y ~art © and c ~gr 6 Xai a. Using our implementation of 
the PF algorithm [5], we can compute a Grobner system for (F’) as follows 


({o-1,a—1],[], [x* + y]) 

([b — 2,4 — 2], [ J, [y, x7}) 

([b — al, [b 2,b—1], [by — y, x7) 

([b — 2}, [a — 2], laxy — 2xy + y, x, yl) 
([],[a—6,b— 2], [axy — bry + by — y, —by + 2? + 2y, 


(2b? — a?b + 2ab? — b? + 2a? — 4ab)y? + (b? — 2b+ 1)y]). 


The above Grobner system has five branches while assuming a = 3 and 
b=2 randomly, the fourth branch corresponds to these values of parameters, 
and so {3xy — 2xy + y,27,y} is a Grébner basis for the ideal (F) |a=3,b=2- 


The question that may arise is: how one could use Grébner system 
to decide the invertibility or singularity of a parametric matrix A? For 
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this purpose, we shall compute a Grobner system {(N;, W;, G;) ‘4 for the 
corresponding parametric linear ideal (Z). Now, for each i if |G;| = |Z| 
then A is invertible provided that N; and W; are satisfied, and it is singular 
otherwise. We illustrate this technique with a simple example. 


Example 3 Consider the following matrix A with parametric entries from 
Kr, s, ¢]. 
r-1 1 ¢t41 


2 2 —Ss 
r+1 3 -l-s 
which corresponds to the following parametric linear ideal: 
T=((r—l)at+yt (t4+1)z,-—sz+ 2x 4 2y, x(r +1) + 8y + (-1- )z). 


A Groébner system of I w.r.t. the product ordering of z ~lex Y ~<lex © and 
t <lex § ~lex T by using our implementation of PF4 algorithm [5] in Maple 
is as follows: 


([], 2+t,r—2r—l1], [(r—-Detyt(tt+lz, 
Gr— dyn (re— st 24 2)e 2 t)z)), 
([r — 1], 2+ ¢], Qn + (—s — 2t—2)z,yt (t+ 1)z, 
(—2—t)z]), ~ ; 
(eed, p= a= 1, (r= 1a + y — 2, Ge Mua ess Dey), 
([r — 2], s+ 2t + 2], a ty +(t+1)z, 2), 
([s+2t+2,r—2], [2+4], eu -1)e, 2 — Oe), 
(2+%,7— 1), IF 2x + (—s+2)z,y—2]]), 
(2+t,s—2,r—2], [], pay all 


This follows the following rank system of A (introduced in [4]): 


iP 2+t,r—2,r—1], 3), 
r— 1, 244], 3), 
2%), r—2,r—1], 2), 
r—Q, 3+2t+2], 2) 
s+2t+2,r—2], [244], 2) 
2+t,r—1], IP 2): 
2+t,s—2,r—2], [], 1). 


? 


o) 


La NN SoS oS 


The rank of a matriz and its determinant are related in the sense that the 
determinant of a square matrix is non-zero if and only if the matrix has full 
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rank. In other words, if the determinant of a square matrix is non-zero, then 
the matrix is invertible and has full rank. Conversely, if the determinant is 
zero, then the matrix is singular and does not have full rank. According to 
this, we can conclude the following system for deciding the invertibility or 
singularity of A. 


([ J, 2+t,r—2,r—1], Invertible), 
([r — 1], 2+], Invertible), 
([2 + ¢], r—2,r—1], Singular), 
([r — 2], s+ 2t+ 2], Singular), 
(Is +2t+2,r—2], [2+¢], Singular), 
((2+¢,r — 1], l, Singular), 
(2+¢,s—2,r—2], [|], Singular). 
For example, if s = —1,r = 2, andt = —3 the fourth branch corresponds 


to these values of parameters, and so Al, , ,-94—_3 18 singular. 
aL eo Sy He} 


3 Inverse Matrix Systems 


The computation of the inverse matrix system via Groébner system has two 
disadvantages: first, when the third component of any triple is “Invertible”, 
the inverse matrix is not calculated and only the invertibility is reported. On 
the other hand, Suzuki and Sato [12] mentioned that the performance of the 
computation of Grébner systems of parametric linear ideals “are generally 
slow and inefficient” (It is worth noting that a parametric linear ideal is 
generated by linear parametric polynomials. In the continuation of this 
paper, a linear parametric polynomial refers to a polynomial that is linear 
with respect to the main variables x, and its parametric coefficients may 
be linear or non-linear). That is why, we apply the efficient algorithm GJS 
to compute the Gauss-Jordan system (without using Grébner systems) for 
parametric matrices based on LDS and GES algorithms proposed in [4]. In 
this direction, we first restate the LDS algorithm to determine the dependency 
of a linear parametric polynomial on a given set of parametric polynomials 
without utilizing the Grébner system and then recall the GES algorithm 
for computing a Gaussian elimination form for a parametric matrix. LDS 
and GES algorithms are crucial in the design of the GJS algorithm for 
computing the Gauss-Jordan form of a matrix with parametric entries. The 
LDS algorithm is essential for solving the problem of parametric linearly 
dependency check. 
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More precisely, let us consider g € S as a parametric linear polynomial 
and a triple (NV, W, G) such that G C S is a linear reduced Grobner basis w.r.t. 
a given monomial ordering according to the constraint sets N,W Cc Ka]. 
The question that may arise is: how to specify the dependency of g on the 
triple (N,W,G)? To answer this question, we described the following LDS 
algorithm in [4]. Below, we use the NORMALForRM function which takes a 
polynomial p, a Grébner basis G = {g1,...,9m}, and a monomial ordering 
~< on the parameters or variables. It returns f and Q = [q,...,qm] such 
that p = qigi +--:+4m9m + f, where f is the normal form of p by G. 


Algorithm 1 LDS (Linearly Dependency System) 

Require: G C S; A linear reduced Grobner basis w.r.t. the product of the 
monomial orderings <x and <a provided that a specification pair (IV, W) 
is satisfied and g € S; a linear polynomial with parametric coefficients 

Ensure: A linearly dependency system of g on (NV, W,G) 

Sys:= {} 
f,Q :=NORMALFORM(g,GROBNERBASIS(N, ~a), <a) 
f’,Q’ :=NORMALFOoRM(f, G, <x) 


if f’ =0 then 
Sys:=Sys UL(N, W, [true, a 0}) } 
else 
A := {a;,,...,@4,} where f! = aj,01, +--+ + 04,24, 


with aj, #0 and xj, >x +++ +x Vi; 
for 7 from 1 tot do 
if a;, is not constant then 
Sys:=Sys 
; U{(N U {ai,, pes Gi; ,},W U {ai;}; [false, Q’, F'lai, =0,...,:,_,=0])} 
else 
Sys:=Sys 
UN U {ai,, es Gis sf; W, [false, Q, F'lai,=0,...,a1,_,=0])} 
Return(Sys) 
end if 
end for 
Sys:=Sys U{(N UA, W, [true, Q’, 0])} 
end if 
Return(Sys) 


The linear dependency system of the following example is computed 
utilizing our Maple implementation of the above algorithm. 
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Example 4 Consider (N,W) = ([m,nJ],|a — 1]) as a specification, G = 
[utv, by—cz, br—acz] a linear Grébner basis, and g = av—bw+mz+nxr—u—y 
a parametric polynomial. We fix the monomial orderings n ~<jex ™ ~lex 
d lex © X<lex b ~<ler @ and w <ler VU ~lex U ~lex 2 Xlex Y ~<lex T ON the 
parameters and the variables, respectively. So a linear dependency system 
of g on G according to (N,W) is as follow: 


({m, n], [c, a — 1], [fatse, |—1, = ,0],-—F + (a+ lv — bw), 
Soe ({m,n,—c],[a—1,a+1], [false,[—1, = 0], &S+(a+1)u— bu)), 

([m,n,—c,a + 1], [O], [ false, [-1, = 0], = — bul), 

([m,n,—c,a+1,0],[ ], [true, [—1 =, 0], 0]). 


The behavior of the above algorithm is illustrated step by step for a sim- 
ple example in [4, Example 15]. Using this algorithm, we re-explain the 
Gaussian elimination for parametric matrices (so-called Gaussian elimination 
systems [4, Definition 10]) and also the efficient GES algorithm to compute 
these systems [4, Algorithm 3]. 


Definition 2 Let A be a parametric matrix over K[a] = K[a1,...,a@m]. A 
Gaussian elimination system of A is a finite set of triples {(.Ni, Wi, Mi), 
such that for any specialization o : K[a] > K, there exists ani satisfying the 
following conditions: 

e The matrix o(M;) is a Gaussian elimination form of o(A), 


e o(p) =0 for each p € N; and o(q) £0 for each q € Wi, 
© Ui VOM) \ Vwew, #) = K™ 


In the following, we restate the efficient GES algorithm for computing the 
Gaussian elimination systems [4, Algorithm 3}. 


Algorithm 2 GES (Gaussian Elimination System) 
Require: M; A matrix with parametric entries over K{a] 
Ensure: A Gaussian elimination system for 
Sys:= {} 
F := the linear system in S corresponding to M 
={(1. 110, FP} 
while A # {} do 
a:= A[l] and A:= A \ {a} 
if a[5] = {} then 
Sys:=Sys U{(a[1], a[2], a[3])} 
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else 
G := al5] \ {a[4]} and g := G/1] and P :=LDS(a[1], a[2], a[3], a[4]) 
for i from 1 to |P| do 
Pli] = (N, W, [flag, Q, f]) 
if flag =true then 
A:= AU{(N, W, a3], g, G)} 
else 
A:= AU{(N,W, a3] U {f},9,@)} 
end if 
end for 
end if 
end while 
Return(Sys) 


Example 5 Let us consider the parametric matrix A in example 3. We get 
the following Gaussian elimination system of A using our Maple implemen- 
tation of the above GES algorithm. 
r—1 1 t+1 
0 2r—4 pense ) 


((], 2+t,r—2,r—1], r—L 

0 0 —t—2 

2 0 -s—2t-—2 
([r — 1], 2+], 0 1 t+1 ), 

0 O —t—2 

r—1 1 —1 
((2+ d], r—2,r—1), | 0 2257 reas? ) 

0 0 0 

1 1 é#+1 
((r — 2], s+ 2t+4 2], c Q =k ) 

0 O 0 

1 1 ¢41 
([s+2t+2,r—2], [2+4], 0 0 -t-2 ]), 

0 0 0 

2 0 -s+2 
((2+t,r— 1], 1, 0 ol —1 ), 

0 O 0 

1 1 -1 
([2+t,s—2,r—2], [], 0 0 O }f). 

0 0 O 
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A matrix is said to be in reduced row echelon form (rref) when it is in row 
echelon form and vectors having one entry equal to 1 and all the other entries 
equal to 0. More precisely, a matrix is in reduced row echelon form if it 
satisfies the following conditions: 


e The left-most nonzero entry in each row (known as the pivot or leading 
entry) is 1, and the column containing this pivot has all other entries 
equal to 0. 


e The pivots in the next rows are to the right of the pivot in the row 
just above. 


e Any zero-row (a row consisting entirely of zeros) is at the bottom. 


The standard algorithm used to transform a matrix into a reduced row 
echelon form is called Gauss-Jordan elimination. However, the Gauss-Jordan 
elimination of a parametric matrix has not, to our knowledge, been studied in 
the literature. On the other hand, many problems in science and engineering 
can be modeled by parametric matrices, and have to be repeatedly solved 
for different values of parameters. So, the Gauss-Jordan of a parametric 
matrix may have a wide range of applications. 

The following example shows that the traditional approach for comput- 
ing the Gauss-Jordan form of matrices may not be used for such a matrix. 


Example 6 Let 


a b+e 0 —1 


where a,b and c are parameters in the real numbers. The reduced row echelon 
form of A computed by the function MTM:-rref(A) of Maple is equal to 


bc—2b—c+3 
1 0 0 a—5+2c 
B= 010 abc—2ab—ac+4a+2c—5 


(a—5+2c)(b+c) 
b-a—b-1 
0 0 1 eres se | 
which is a matrix of rank 3 for all values of parameters. However, this is 


undefined for the values a = 3,b = 2 andc = 1, where the Gauss-Jordan 
form of A|,a=3,=2,c=1} 18 the following Matrix with rank 2. 
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In fact, one may be interested in finding for what values of parameters, the 
Gauss-Jordan form of the above matrix is of rank 2 and for what values it is 
of rank 3. In doing so, we introduce the notion of the Gauss-Jordan system 
for a parametric matrix, and we present an algorithm for computing it. 


Definition 3 Let A be a parametric matrix over Klay,...,@m]|. A triples 
set {(Nj, Wi, M;)}{_, is called a Gauss-Jordan system of A if 


¢ UV) \ VIuew, #) =K 


e for each specialization o : K[a| > K, there is ani such that o(M;) is 
a Gauss-Jordan form of o(A). 


mm 


Using the LDS and GES algorithms proposed in [4], we present an efficient 
algorithm to compute a Gauss-Jordan system for parametric matrices. Below, 
given an n X m parametric matrix M over Klal], we associate a linear 
parametric polynomial system FC S. 


Algorithm 3 GJS (Gauss-Jordan System) 
Require: M; a parametric matrix over K[a] 
Ensure: A Gauss-Jordan system for M 
Compute a Gaussian Elimination System {(N;, Wi, M;)}{_, 
for M using GES algorithm [4] 
Sys:= { } 
for 7 from 1 to @ do 
B; := the reduced row echelon form of M; 
Sys :=Sys U {(Ni, Wi, Bi)} 
end for 
Return(Sys) 


Theorem 4 The GJS algorithm terminates and accurately computes a 
Gauss-Jordan system. 
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Proof: The termination and correctness of the GES algorithm ensure 
the termination and correctness of the GJS algorithm. Specifically, for each 
triple (N;,W;, ;), the first entry of any nonzero row from the left of M; is 
nonzero according to specification (N;,W;). Thus, the Gauss-Jordan form 
can be completed for each triple. Of course, if a row has all zero entries, then 
it doesn’t have a leading entry. Therefore, Sys is a Gauss-Jordan system 
for M at the end. 


Example 7 Consider the matrix A mentioned in example 3. Using our 
Maple implementation of the GJS algorithm we obtain the following Gauss- 
Jordan system of A. 


1 0 0 
(LI, Dart 21) 0 1 0 ); 
00 1 
1 0 0 
(ree; [2 + ¢], i ); 
00 1 
—2 
1 0 ez 
—s-2 
([2 + t], [(r —2,r— 1], OS A eee Is 
0 0 0 
1 1 0 
({r — 2], [s + 2é + 2], 00 14), 
0 0 0 
1 1 0 
(Is +2t+2,r—2], ([2+¢], OO Ne 
0 0 0 
10 3841 
(249 =1], I Os a = )s 
| 0 0 0 | 
11 -l 
(2s = 2er= 2), [], 0 0 0 ) 
0 0 O 
For example, the parameter values r = 3,t = —2, and s = —4 hold at the 


parametric constraints of the third branch of this system, and so the reduced 
row echelon form of A|p=3t=—2,s=—4 18 as follow: 
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=2 

2 A Lt) 33 iG. cate 

~s—2 

Emel 20 22 2i4ey| |e | oe es ives Oe es Soe 
4 3 3 0 0 O 00 0 


r=3,t=—2,s=—A4 


In the following, first, we will briefly provide an overview of the Gauss-Jordan 
method for computing the inverse matrix with constant entries. We then 
introduce and compute an inverse matrix system utilizing the extension of 
the Gauss-Jordan method (the Gauss-Jordan systems) to handle matrices 
with parametric entries. 

To find the inverse of a n x n matrix A: adjoin the identity matrix I to 
the right side of A, thereby producing a matrix of the form [A|/] and then 
apply the row operations to this matrix until the left side is reduced to J. If 
successful, these operations will convert the right side to A~', so that the 
final matrix will have the form [J|A~1). 

In many scenarios, the elements of a matrix are not fixed constants but 
instead depend on parameters or variables. These matrices with parametric 
entries play a pivotal role in modeling systems subject to change or vari- 
ability. The ability to compute the inverse of such matrices is essential for 
solving systems of linear equations and gaining insights into the behavior 
of parametric systems. So, applying the Gauss-Jordan method to compute 
the inverse of matrices with parametric entries is an appropriate idea. In 
this direction, we introduce the concept of an inverse matrix system for 
parametric matrices. Then, we present the IMS algorithm, which extends 
the Gauss-Jordan method, to compute such systems. 


Definition 5 Let A be a square parametric matrix over Kia] = K[a1,...,@m| 
and A =det(A). A finite triples set {(N;, Wi, Ai)}$§_, U{(N’;, Wi Ai Veen 
is an inverse matrix system of A if for each specialization o : Kia] > K 


there is ani € {1,...,0} or one 7 € {€+1,..., 0} such that 
e o(A;) is inverse matrix of o(A) (or Al, = Singular) 


e o(p) =0; Vp € N; C K[a] (or Vp € Nj Cc K[al) 


o(q) #0; Vqe W; C Kfa] (or Vg € W; Cc Rial) 


Ure VIN) \ VT wew, &) = K™ \ V(A) 
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For computing the inverse of a parametric matrix Anyn, we let Sys be a 
variable which is initialized to an empty set, and finally, it is the output 
inverse matrix system. Similar to matrices with constant arrays, at the 
beginning we can adjoin the identity matrix In.» to the right side of A, 
forming the matrix [A|/] and then apply the GES algorithm to this matrix 
for computing a Gaussian elimination system of [A|I]. There are two kinds 
of triples in this step: those in which the rank of the produced matrices is 
less than n when the state is “Singular”, and those whose matrices have 
rank n, which are the candidates of the “Invertible” state. The goal is to 
compute inverses, so the Gauss-Jordan systems are computed only when 
the rank equals n. Therefore, for the first kind, without any additional 
computation, the triple (N;,W;, “Singular” ) is added to Variable Sys and 
the GJS algorithm idea only is applied to the matrices of the second kind 
of triples to compute a Gauss-Jordan system of [A|J]. At the end of this 
procedure, for the second triples kind, the right side of the produced matrices 
n x 2n is A~!. Briefly, each recorded triple in Sys is of the form (N, W, 
Singular or A~') where (N,W) is a conditions pair and A~! is the inverse 
matrix of A according to (N, W) in non-singular cases. 


Algorithm 4 IMS (Inverse Matrix System) 


Require: A,.,; a parametric matrix over K{a] 
Ensure: An inverse matrix system for A 
ME on = [Anxn|Inxn] 
{(Ni, Wi, M;)}{_, := A Gaussian Elimination System 
for M using GES algorithm 
Sys:= { } 
for i from 1 to ¢ do 
if rank(M;) <n then 
Sys :=Sys U {(M, Wi, “Singular” )} 
else 
B; := the reduced row echelon form of M; 
R; := The right sub-matrix of B; with order n 
Sys :=Sys U {(Mi, Wi, Ri)} 
end if 
end for 
Return(Sys) 
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Theorem 6 IMS algorithm terminates after a finite number of steps and is 
correct. 


Proof: The termination of the algorithm is ensured by the finiteness of 
the GES algorithm. The correctness of the algorithm is guaranteed by the 
correctness of the LDS algorithm and the authenticity of the Gauss-Jordan 
method for the computation of inverse matrix with constant arrays. Indeed, 
by the structure of the algorithm, we explore applying the Gauss-Jordan 
method to compute the inverse of matrices with parametric entries. More 
precisely, it discusses parameter space using the LDS algorithm and for all 
consistent possible states decides the invertibility or singularity of matrices 
according to the parametric constraints. If the left sub-matrix M is equal to 
the identity matrix J then the triple (N;,W;, R;) is added into the global 
variable Sys where R; is the right sub-matrix M. Otherwise, the triple 
(N;, W;, “Singular” ) is added into the global variable Sys which means that 
A is singular with respect to the corresponding conditions pair (N;, W;). 


We illustrate the behavior of this algorithm with a simple example. 


Example 8 Let us consider the parametric matrix from example 3 in the 
following. 


r—-1l 1 t+4+1 
2 2 —s 


a 
I 


r+1 3 -l-s 


In the initial stage of the IMS algorithm, a Gaussian elimination system of 
A is computed (see example 8.3). The rank of the corresponding matrices 
of the last five triples of this system is less than 3. As a result, without 
any additional computation, the triple {(N;, Wi, “Singular” )}?_, is added to 
Variable Sys. Then, the GJS algorithm idea is applied only to the first two 
triples for completing the Gauss-Jordan system for [A|I]. Finally, using our 
Maple implementation of the IMS algorithm an inverse matrix system of A 
is obtained as follows: 
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‘Parametric Constraints 


Inverse Matrix System (IMS) 


W,=([] 


s—2 s+3t+4 s+2t+2 
1/2 (€-£2)(r—2) 1/2 CaN) 1/2 GG 
M > [ 1, rs—s—2 rstrt+2r—s+t rs—s+2t+2 
Wi =(24+t,r—-2,r—-] 12 Gey V2 Gaya) 1/2 G3) 
eo) Gage pe tee ae 
1/2 a 1/2 otsits 1/2 sie 
No = [r— 1, -1 t+1 t41 
We = [24+ 1] a ~ t+2 42 
Co) eee (eo) ee (oe) a 
Nz = [2+ ¢], : 
Wea aoe Singular 
Na = [r ~~ 2], : 
Wa = [s+ 2t+2] Singular 
Na le penal, Singular 
Ws =([2+¢] 
i Z Weed, Singular 
Nam Oat ene .e) Singular 


The computation of the inverse of parametric matrices (inverse matrix 
systems) can be used to solve systems of linear equations involving parametric 
coefficients. This opens up new possibilities for analyzing and predicting 
the behavior of systems with varying parameters. There are many practical 
examples from engineering, physics, and economics where matrices with 
parametric entries and their inverses play a critical role. These examples 
showcase the versatility and applicability of the Gauss-Jordan method in 
parametric modeling relevance in addressing real-world challenges. 
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4 Algorithm Performance 


In this section, we provide the performance of computing the inverse of 
parametric matrices using the Gauss-Jordan method, including the used 
memory and the timing duration of the computations via some examples. 
For this purpose, we have implemented all the algorithms described in this 
paper in Maple 18. The following parametric matrices are chosen from the 
ring Kja, b,c,d,m,n], and our aim is to compute an inverse matrix system 
of each matrix. 

The results are shown in the tables below, where the timings were 
executed on a personal computer with Ryzen 6800, 8 GB RAM, and 64 
bits under the Windows 10 operating system. The third and fourth columns 
in these tables indicate the CPU time (in seconds) and the amount of 
required memory (in gigabytes) of the procedure. The last column shows 
the number of triples of the computed IMS algorithm. To illustrate the 
efficiency of our algorithm, we compare the performance of two methods, 
i.e. Grébner systems and the IMS algorithm. However, this comparison is 
somewhat unfair because the IMS algorithm not only reports invertibility 
but also computes the inverse matrix in that case. In this direction, suppose 
that the following F,,...,F7 are the parametric linear polynomial systems 
corresponding to the above matrices in the parametric polynomial ring 
K[a, b,c, d,m, n][#1, 22, £3, 24, U5, LG, L7, Lal. 

e Fi = [(a—1)a, +22 4+ (4+ 2)az, (b+ 1)xe — 23,341 + (a — b)xy — 223] 
e Fy = [(a— 1), +242 + (a4 Lag — v4, 2a, + (b+ 1)x3 — 44,2104 (C+ b)to + 
(a — b)ag3 — 2a4, (a — 1)a1 + (b— c)x2 + (a — c)z3] 


e F3 = ((2—m)a, + nx 4+ (c— 1)a3 — @4, —ba4 — 22 + %3,MxX, + CHQ + (1 -— 
b)a3 + 2x4, —X1 — cXg + (a — 1)2°3] 


e Fy = [nx + x3 — U5, bay + (a — L)ae + og + 205, —201 + (m— 1)ag +43 + 
2x5, (nm — 1l)ag + 3x3 + bag, —Xq + (14 d)axzg + 224] 


e Fs = ((m—1)a, +43 + 245, nx +243 +04 — 46, —221 + (a— 1)x3 — 25, dxr4— 
t+ 2X5 — 6, —X3 + (b— 1)u5, carp + 301 + VQ — 25] 

e Fe = [zo +(m—1)a3+(m+1)a5, (m? —1)a1 + 2x3 - v4 +264+M2x7,22+(n—- 
1)zg + (n+ 1)a4, —ro + (m + 2)r3 — 26, -41 + (—M +4 1)ze — 24 +25, —-22 + 
223 + (m — 2)a7, (n — 2)23 +24+(m—n)as + (m? — 4) a6) 


e Fy = [(a—1)a, +23 4+ 204+(at1)25 4+ (a? +a)r7, 22+ (a—2)x3 +(a+3)a5 + 
x7, (a? —1)a, + 223 — 244+ 4504-26 +28, %2 + (a—1)x3 + (14+ 3a)x4 +4 (a* 
1)zg, (a* — a*)x1 — xo + 23 + 2arg — 2g,—21 + (4 + 1) xq — 24 4+ 25, —22 + 
273 + x5a + (a — 1)a7, -2r3 + v4 + (a? + 1)a5 + (a? — 1)x6 + (a® — a?) 27] 
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Time Used 
Matrix (Sec) Memory | Branch 
eC (GB) 
1 c+2 
+1 -1 0.23 0.02 6 
a—b —2 3x3 
a+1-1 
0 b4+1-1 
0.98 0.14 15 
c+ba—b —2 
a—-1lb-—-ca-c 0 AXA 
aeons n c-1-l1 
1 —-b 
2.01 0.34 21 
c 1-b 2 
—-ca-1 0 AXA 
ro 0 on 1 oes 
b a-1l 0 12 
—2m—-1 1 0 2 4.67 0.79 45 
Oo n-1 3 +060 
0 l+d20 4]... 
oxo 
m—-1 0 1. 0 
Te Oe 21 
-1 0a-10 -1 0 
6.78 1.02 66 
Oo di -tl 
-1 0b-1 0 
0 0-1 ie 8x6 
| 0) 1 m-1 0 m+4+i1 0 0 ] 
2-1 0 ) -1 0 1 m 
0 1 n-lil+n 0 0 0) 
0 -1 m+2 0 0 —1 0 6.28 0.97 28 
-l1 1l-m 0 -1 1 0 0 
0 -1 2 0 0 0 m—2 
2 
0 0 n—-2 1 m-nm*-4 0 7x7 
a—1l 0 1 2 atl 0 a@+a 0 
0 1 a-2 0 a+3 0 1 0 
a®—1 0 2 -1 a 1 0 1 
0 1 a-11+3a 0 0 0 at-1 
ol 0.45 8 
a*—a? -1 2 2a 0 HL 0 0 
-1 a+l1 0 -1 1 0 0 0) 
0 -1 2 0 a 0 a-1 0 
0 0 -2 1 es 0 = OP Re AU ee erat 0 decal) tte = | 0 
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The efficient PGBMAIN algorithm can now be used to compute Grobner 
systems. However, a drawback of this algorithm is that the output Grobner 
system may contain several branches, resulting in the corresponding Grobner 
basis being {1}. By consolidating these branches into a single branch, the 
consistency check, as well as the used memory and time of computations, 
can be reduced, significantly improving the algorithm’s performance. This 
simple modification was included in the PGBMAIN algorithm in [6] and we 
use it to compute Grodbner systems of the seven lists mentioned above. 


Example | Time (sec.) | Used Memory (GB) 
Fy 0.28 0.02 
Fo 1.12 0.17 
F3 3.64 0.63 
Fy 3.93 0.56 
Fs 6.56 1.28 
F¢ >300 — 
F7 29.65 5.26 
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